Existence of One-Body Barrier Revealed in Deep Sub-Barrier Fusion 
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Based on the adiabatic picture for heavy-ion reactions, in which the neck formation in the one-body system 
is taken into account, we propose a two-step model for fusion cross sections at deep subbarrier energies. This 
model consists of the capture process in the two-body potential pocket, which is followed by the penetration of 
the adiabatic one-body potential to reach a compound state after the touching configuration. We describe the 
former process with the coupled-channels framework, while the latter with the WKB approximation by taking 
into account the coordinate dependent inertia mass. The effect of the one-body barrier is important at incident 
energies below the potential energy at the touching configuration. We show that this model well accounts for the 
steep fall-off phenomenon of fusion cross sections at deep subbarrier energies for the 64 Ni+ 64 Ni and 58 Ni+ 58 Ni 
reactions. 

PACS numbers: 25.60.Pj, 24.10.Eq, 25.70Jj 



Heavy-ion fusion reactions at low incident energies provide 
a good opportunity to study the quantum tunneling phenom- 
ena of many-particle systems. Because of a strong cancel- 
lation between the repulsive Coulomb interaction and an at- 
tractive short range nuclear interaction between the colliding 
nuclei, a potential barrier, referred to as the Coulomb barrier, 
is formed, which has to be overcome in order for fusion to 
take place. In heavy-ion reactions, because of a strong absorp- 
tion inside the Coulomb barrier, it has been usually assumed 
that the compound nucleus is automatically formed once the 
Coulomb barrier has been overcome. The coupled-channels 
(CC) approach based on this picture has been successful at 
energies close to the Coulomb barrier, where the inner turning 
point of the Coulomb barrier is well outside the touching point 
of the colliding nuclei fljj] . 

Recently, fusion cross sections have been measured for the 
first time at deep subbarrier energies for medium-heavy mass 
systems, such as 64 Ni+ 64 Ni, 58 Ni+ 58 Ni and 64 Ni+ 89 Y Hi]. 
The experimental data indicate that fusion cross sections fall 
off much faster than the exponential energy dependence ex- 
pected from a usual tunneling picture, as the incident energy 
decreases. Although it has been argued that this hindrance of 
fusion cross sections may be explained if one phenomenolog- 
ically introduces a considerably diffuse nuclear potential yj, 
the physical origin of the steep fall-off phenomenon has not 
yet been understood (see also Ref. [i]). 

At energies well below the Coulomb barrier, the inner turn- 
ing point is comparable to, or even smaller than, the touch- 
ing point. In that situation, the frozen density approximation, 
which has often been employed in constructing the internu- 
cleus potential [6], breaks down, and one has to treat explicitly 
the dynamics after the touching configuration. In this connec- 
tion, Misjcu and Esbensen have recently proposed a potential 
energy with a shallow pocket based still on the frozen den- 
sity approximation That is, the outer region of the po- 
tential is constructed with the double folding procedure [6], 
while the phenomenological repulsive core due to the satu- 
ration property of nuclear matter is taken into account in the 
inner region 0,1a] • It was shown that the CC calculation with 



such shallow potential well reproduces the steep fall-off phe- 
nomenon for the 64 Ni+ 64 Ni reaction 



The approach of Migicu and Esbensen is based on the sud- 
den picture for nuclear reaction, that is, the reaction takes 
place so rapidly that the colliding nuclei overlap with each 
other without changing their density. However, it is not obvi- 
ous whether the fusion dynamics at deep subbarrier energy is 
close to the sudden limit or to the adiabatic limit, where the 
nuclear reaction is assumed to take place much more slowly 
than the dynamical density variation of colliding nuclei. Since 
one would not know a priori which approach is more reason- 
able, it is important to investigate both the possibilities |9(]. 

In this paper, we investigate the adiabatic approach in ex- 
plaining the steep fall-off phenomenon of fusion cross sec- 
tions. Notice that both the sudden and the adiabatic ap- 
proaches would lead to a similar result to each other in the 
region where the colliding nuclei do not significantly over- 
lap. Our model here is to consider the fission-like adiabatic 
potential energy surface with the neck configuration after the 
colliding nuclei touch to each other. This one-body potential 
acts like an inner barrier which has to be overcome to reach 
the compound state. It is this residual effect which we would 
like to discuss in connection to fusion cross sections at deep 
subbarrier energies. 

In order to illustrate how the adiabatic approach works, 
Fig. [TJ shows the potential energy for the 64 Ni+ 64 Ni reaction 
obtained with the Krappe-Nix-Sierk (KNS) model llioll as a 
function of the center-of-mass distance R. In the KNS model, 
the saturation property of nuclear matter is phenomenologi- 
cally taken into account. It has also been shown that the KNS 
model is consistent with the potential obtained with the energy 
density formalism with the Skyrme SkM* interaction fiUl . 
The parameters in the KNS model are taken as ao=0.68 fm, 
a. v =21.33 MeV and K s =2.378 from FRLDM2002 |Qj]. The ra- 
dius parameter is fine-tuned as ro = 1 .204 fm in order to fit the 
experimental fusion cross section at high incident energies. 
The touching configuration is denoted by the filled circle in 
the figure. For distances larger than the touching point, the po- 
tential energy for the two-body system is calculated as the sum 
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of the Coulomb energy for two point charges and the nuclear 
energy given by Eq. (17) in Ref. 11 Oil . For the one-body sys- 
tem after touching two nuclei, we assume that the shape con- 
figuration is described by the Lemniscatoids parametrization 
(see the inset in the figure) fl3ll . and calculate the Coulomb 
and surface integrals for each configuration ifioll . 

We find that the value of the potential energy at the touching 
configuration V t ouch is 88.61 MeV. This is exactly the energy 
E s at which the experimental fusion cross section start to fall 
off abruptly in this reaction yj]. This strongly suggests a cor- 
relation between the observed fusion hindrance and a process 
after the two nuclei overlap each other. For a comparison, the 
sudden potential which Migicu and Esbensen considered |7|] is 
denoted by the dotted line in the figure. We find that the adia- 
batic KNS potential and the sudden potential almost coincide 
with each other outside the touching radius. 

In order to describe the two-body process from a large dis- 
tance to the touching point, we employ the standard CC for- 
malism by taking into account inelastic excitations in the col- 
liding nuclei. However, it is not straightforward to extend this 
treatment to the one-body process. In the CC formalism, the 
total wave function is expanded with the asymptotic intrinsic 
states of the isolated nuclei, in which one usually restricts the 
model space only to those states which are coupled strongly 
to the ground state. Apparently, such asymptotic basis is not 
efficient to represent the total wave function for the one-body 
di-nuclear system, and in principle one would require to in- 
clude all the intrinsic states in the complete set. This is almost 
impossible in practice. Moreover, the adiabatic one-body po- 
tential with the neck configuration already includes a large 
part of the channel coupling effects, and the application of the 
standard CC formalism would result in the double counting of 
the CC effect. 

In order to avoid these difficulties, we here propose a sim- 
ple phenomenological model, in which the two- and one-body 
processes are defined independently and time-sequentially. 
The fusion cross section in this two-step model then reads 

7ih 2 

<j(E) = —£(2£+l)T e (E)P lM (E,£), (1) 

where ji and E denote the reduced mass and the incident en- 
ergy in the center-of-mass system, respectively. 7) is the cap- 
ture probability for the two-body process estimated with the 
CC method. Pibd is the penetrability for the adiabatic one- 
body potential to reach the compound state after the touching 
of two-body potential, which plays an important role at ener- 
gies below Vtouch (i-e., below the dashed line in Fig. [TJ. At 
these energies, the fusion reaction is not described only by 
the two-body potential, but the potential which governs the 
fusion dynamics is switched from the two-body to the adia- 
batic one-body potential at the touching configuration. Only 
after overcoming (or penetrate through) these two- and one- 
body barriers, the system can form a compound nucleus. One 
may regard the one-body penetrability Pibd as a fusion spec- 
troscopic factor, which describes the overlap of wave function 
between the scattering and the compound states. 

In order to estimate the capture probability T( within the 
two-step model, we cut the two-body potential at the touch- 
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FIG. 1: (Color online) One- and two-body potential energies for 
64 Ni+ 64 Ni obtained with the KNS model as a function of the center- 
of-mass distance. The shape for the one-body configuration de- 
scribed by the Lemniscatoids parametrization is also shown. The 
filled circle and square denote the touching configuration and the 
ground state of the compound system, respectively. The dotted lines 
is the sudden potential taken from Ref. [7|]. 



ing configuration as shown in the upper panel of Fig. [2] The 
capture probability does not depend strongly on how to cut 
the potential, since only the lowest two-body eigen potential, 
which is obtained by diagonalising the coupling Hamiltonian 
id 0, Eltl, is relevant at deep subbarrier energies. As indi- 
cated by the dashed line in the figure, the inner turning point 
for the lowest eigen potential is still far outside the touching 
distance. Thus, the actual shape of the original potential in 
the inner-barrier region influences little on the penetrability. 
Another view is that the incoming wave boundary condition 
(IWBC) is imposed in the CC calculation at the touching dis- 
tance so that the capture probability is defined at the touching 
configuration, although in the actual calculations we impose 
the IWBC at a distance somewhat smaller than the touching 
point in order to avoid the numerical error. For simplicity, we 
employ a sharp cut-off of the two-body potential in this paper. 

In order to estimate the one-body probability Pibd, we use 
the WKB approximation. We assume that the reflected flux in 
this process does not return to the two-body system, but exits 
through the muti-dimensional potential energy surface in the 
one-body system. The penetrability then reads P\bd(E,£) = 
e -2S{E,e)^ yyjjgj-g S(E,£) is the action integral with the coordi- 
nate dependent inertia mass M(R), 



S(E,£) = J R b dR^?^(E-V lbd (R,m- (2) 

Here, R a and R^ are the inner and the outer turning points, re- 
spectively (see the lower panel of Fig.|2|. V^d is the adiabatic 
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FIG. 2: (Color online) The intemucleus potential used in the two-step 
model. The solid line in the upper panel denotes the KNS potential 
for the two-body process, which is cut at the touching configuration, 
while the dashed line denotes the lowest two-body eigen potential. 
The dash-dotted line denotes the position at which the in-coming 
wave boundary condition (IWBC) is imposed in the CC calculation. 
The solid line in the lower panel denotes the adiabatic one-body po- 
tential inside the touching distance. 




FIG. 3: (Color online) Fusion cross sections for the Ni-H Ni reac- 
tion calculated with the two-step model. The filled circles denote the 
experimental fusion cross section, taken from Ref. |2]. The solid and 
dotted lines denote the fusion cross section obtained with the normal- 
ization factor for the mass inertia of k = 46 and 0, respectively. The 
dashed line denotes the corresponding capture cross section. The 
dash-dotted line is obtained with the Woods-Saxon potential, while 
the dash-dot-doted line shows the result in the absence of the channel 
coupling effect. 



one-body potential energy given by 

V lhd (R,i) = V C (R)+V S (R) + e{e + { ^ 2 +^E R , (3) 

where Vc, Vs and / are the Coulomb and the surface energies 
and the moment of inertia for the rigid body, respectively. Er 
denotes the centrifugal energy at the touching configuration. 
Note that the last term in Eq. (O comes from the conservation 
of the energy and angular momentum between the two- and 
one-body systems in the sticking limit llrjll . 

We now apply the present two-step model to the fusion re- 
action of the 64 Ni+ 64 Ni system. To this end, we use the KNS 
potential energy already shown in Fig. Q] In the energy re- 
gion discussed in this paper, we expect that the the Lemnis- 
catoids parametrization provides a reasonable approximation, 
because the neck formation is still small as shown in the in- 
set of Fig. Q] This parametrization has an advantage in that 
the configuration is described with only one parameter for a 
symmetric system. In addition, one obtains a smooth connec- 
tion between the one- and two-body potential energy curves, 
since the change of the configuration shape across the touch- 
ing point is rather natural. As for the inertia mass M, we take 
the linear combination between the irrotational-flow mass in 
the Werner- Wheeler approximation [17], My, and the reduced 
mass, pL. That is, M(R) = k (Mq(R) — ju) + ju, where k is the 
normalized factor. The renormalization factor is necessary, 
since the liquid drop model with the irrotational-flow mass 
Mq overestimates the vibrational excitation energy hcoo for the 
first 2 + state fl8il . In the calculations presented below, we use 
the normalization factor, k = 46, which leads to the vibrational 
energy of 0.2 TicOq. Notice that the inertia mass M is in agree- 



ment with the reduced mass ji at the touching configuration. 

In order to compute the capture probability T( with the CC 
framework with a sharp-cut KNS potential, where the form of 
the coupling potential is not known, we modify the computer 
code CCFULL [19] and estimate the nuclear coupling term 
with the numerical derivative of the nuclear potential up to 
the second order. The coupling scheme included in the calcu- 
lations, as well as the deformation parameters, are the same as 
in Ref. [2]. To be more specific, we include the coupling to the 
low-lying 2 + and 3~ phonon states, two-phonon quadrupole 
excitations, and all possible mutual excitations both in target 
and projectile nuclei. The potential depth in the inner-barrier 
region for the sharp-cut KNS potential and the position of the 
IWBC are chosen as V Q = 70 MeV and R mBC = 8.0 fm, re- 
spectively. These values are determined using the the Woods- 
Saxon (WS) potential with V W s=75.98 MeV, rws=l-19 fm, 
and fl\vs=0.676 fm. We have checked the numerical stability 
of the calculations at extremely low incident energies by com- 
paring the obtained result with the one in the multi-channel 
WKB approximation 12011 . 

Figure 3 shows the fusion cross sections thus obtained. It is 
remarkable that the fusion cross section obtained with k — 4-6 
for the coordinate dependent mass is in good agreement with 
the experimental data (see the solid line). The corresponding 
capture cross sections, obtained by setting Pibd = 1 in Eq. (Q]), 
is denoted by the dashed line. As a comparison, the result 
with the WS potential is also shown by the dash-dotted line. 
We see that the discrepancy between the capture cross section 
obtained with the WS potential and the experimental data is 
improved by taking into account the saturation property sim- 
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ulated by the KNS potential, and a further improvement has 
been achieved by taking into account the one-body barrier in- 
side the touching configuration. The result with k = is de- 
noted by the dotted line. The difference between the solid and 
the dotted line is small, indicating the negligible effect of the 
coordinate dependence of mass inertia in the energy region 
discussed in this paper. We have applied the two-step model 
also to the 58 Ni+ 58 Ni system. We found that the agreement 
with the experimental excitation function B21I1 is as good as 
for the 64 Ni+ 64 Ni system shown in Fig. [3] 

The present two-step model is in the opposite limit to the re- 
cent sudden model of Migicu and Esbensen |7|,|8|]. As long as 
the fusion cross sections are concerned, both the models pro- 
vide similar results, at least for the 64 Ni+ 64 Ni reaction. How- 
ever, the origin for the fusion hindrance is different between 
the two approaches. In our two-step model, the fusion hin- 
drance takes place due to the penetration of the inner one-body 
potential. On the other hand, in the sudden model, which uses 
a shallow potential, the hindrance occurs because of the cut- 
off of the high angular-momentum components in the fusion 
cross section. The average angular momentum of the com- 
pound nuclei estimated with the sudden model would there- 
fore be much smaller than that of the present adiabatic model. 
It is thus interesting to measure the average angular momen- 
tum of the compound nucleus at deep subbarrier energies, in 
order to discriminate the two approaches. 

We would next like to comment on the recent experimental 
data for le O+ 197 Au, where the fusion hindrance was not ob- 
served li22ll . We estimate the potential energy at the touching 
configuration, Vtouch, to be 68.23 MeV if we use ro=1.2 fm 
in the KNS potential. This is nearly equal to the lowest in- 
cident energy performed in the experiment. Thus, the fusion 



cross sections have to be measured at lower energies in order 
to observe the fusion hindrance for this system, as has been 
speculated in Ref . i22ll . 

To summarize, we have proposed the adiabatic two-step 
model for fusion cross sections at deep subbarrier energies. 
By applying this model to the 64 Ni+ 64 Ni and 58 Ni+ 58 Ni reac- 
tions, we have shown that the penetration of the adiabatic one- 
body potential with the neck configuration after the touching 
of two colliding nuclei is responsible for the steep fall-off of 
fusion cross sections observed recently in the experimental 
data. The effect of the one-body potential is important only 
at energies below the potential energy at the touching config- 
uration. In this way, the two-step model provides a natural 
origin for the threshold energy of fusion hindrance discussed 
inRefs. JlH. 

In Ref. J3l, it was shown that the experimental fusion cross 
sections for the 16 O+ 208 Pb system follow the exponential en- 
ergy dependence at deep subbarrier energies. This is in con- 
trast to the behaviour in the medium-heavy systems discussed 
in Refs. J3,[3t]. It would be an interesting future work to ap- 
ply the present two-step model to this reaction and to clarify 
the difference between the mass asymmetric and symmetric 
systems. 
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